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ABSTRACT 

This  paper  presents  abstract  error  estimates  for  mixed  methods 
for  the  approximate  solution  of  elliptic  boundary  value  problems. 
These  estimates  are  then  applied  to  obtain  quasi-optimal  error 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  presents  abstract  error  estimates  for  mixed  methods  for  the 
approximate  solution  of  elliptic  boundary  value  problems.  In  a mixed  method, 
one  introduces  an  auxiliary  variable,  usually  representing  another  physically 
important  quantity,  in  order  to  write  the  differential  equation  as  a lower 
order  system.  One  then  considers  Ritz-Galerkin  approximation  schemes  based  on 
a variational  formulation  of  this  lower  order  system.  The  abstract  estimates 
are  then  applied  to  obtain  quasi-optimal  error  estimates  in  the  usual  Sobolev 
norms  for  four  examples:  three  mixed  methods  for  the  biharmonic  problem  and 
a mixed  method  for  second  order  elliptic  problems.  In  the  biharmonic  problem, 
for  example,  one  obtains  by  these  mixed  methods  direct  approximations  to  the 
stream  function  and  vorticity  in  hydrodynamical  problems  or  to  the  displace- 
ments and  forces  (or  moments)  in  elasticity  problems. 


Accession  For 


NT  IS  GRA&I  (D 
DDC  TAB  j* 
Unannounced 

Justification 

r 

By 

Pi  y * 1 f>r  f 

A vp  i 

^ i * 't  C •-"‘‘•s 

Dist 

f\ 

A , all  and/or 

special  , 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ERROR  ESTIMATES  EUR  MIXK1'  METHODS 


K.  S.  Falk  and  J.  F.  Osborn 


1.  Introduction 


In  l SI  llrern  studied  Ri  t z-dal erkin  approximation  of  saddle-point  ptol'Kw  aiising 


in  connection  with  tagiange  multipliers.  These  problems  have  the  tv 


Given  f < v*  and  g t w’,  find  (u,\|)  . V * w satisfying 


j a(u,v)  t h(v,y)  - (f,v)  v v , V 

^ h(u , »■)  - (g , „•)  V y-  , w 


who  it*  V and  W ate  real  Hilbert  spaces,  and  a(*,*)  ami  !>(•,•)  .no  Ustnded 


bil  inoat  forms  on  V ' V and  V * w,  respect  ivel  y . 


Given  finite  dimensional  spaces  V.  * v and  w.  w,  0 * h \ , the  Kit  ? -Galer K tn 

h It 


approximation  ( '*  • W ) is  the  solution  of  the  following  problem  t 


Find  (u  , ) ♦ V " W satisfying 

it  h It  h 


a < \i  , v)  ♦ b(v,4>.)  - (f,v)  V v f V. 
h h It 


b(u  ,v')  - (g,c)  V v*  « Vi 
n it 


The  matoi  assumptions  in  Bter.zi'a  tesultsaie 


"l,l'  ‘\lvl|V  ' \i  11 1,11  V V u ' 2h  v ''  . 

V * Z.  V 


where  > ' 0 is  a constant  independent  of  h ,„,d 


tv  V.  htv.v 

It 


0 v , w^!,  and 


k , Hell 
0 


"w  v C • Wh  and  v h 


where  ' 0 is  independent  of  h . 


lising  (I."  and  (1.4'  Hrerri  pit'ves  the  fol  lt*wing  ettot  estimate  toi  the  ,ppi«'\\m.it 


met  hod  det  outlined  by  (I..*': 


department  of  Mat  hrmat  it's , Nutgets  Pnivetsity,  New  IMtinswioV.  N d.  and 
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j 


1 


(1  s)  llu-u  II  t Hi,-*. II  ■ 0(inf  llu-ill  * inf  <1  tf-- nil  > V h . 

h \ h w v W 


•i<  W, 


In  11..')  KabuSka  studied  Kitr  dalerkin  approx imat  ion  of  general  var iat ionally  posed 
problems.  The  main  result  of  (1,2),  applied  to  (l  .11  and  (1.2),  is  that  (1.5)  holds 
pi  ov ided 


1 1 .el 


nip 


(v,*).v)i*wh 


a (u ,v'  j b(v,j)  * b(u |£ )J 

..  ....  ~ - * 

\ 


B v II  . ii*-il  o(llullv  4 IUV  v ‘“•"'’'W"'1  v h. 


where  i , 


0 is  independent  of  h . 


It  is  oleai  from  (1,2,51  that  (1.0  and  (1.4)  hold  if  and  only  if  (1.6)  holds.  (1.1)- 
(1.4)  or,  equivalent  ly,  (l.e)  is  referred  to  as  the  stability  condition  foi  this  approxi- 
mat ion  method. 

The  results  of  11,2,51  can  he  viewed  as  a strategy  for  analyzing  these  approximation 
methods:  the  approximation  meth«xl  is  characterized  by  certain  bilinear  forms,  norms  (spaces), 

and  families  of  finite  dimensional  approximating  spaces,  and  if  the  method  can  be  shown  to 
be  stable  with  respect  to  the  chosen  norms,  then  the  error  estimates  in  t lies.'  norms  follow 
in  a simple  mantlet  provided  the  bilinear  foims  are  Ixmnded  and  the  approximation  properties 


of  V and  W are  known  m these  nouns.  These  results  can  be  used  t o analyze,  for 
h h 


example,  certain  hybrid  methods  fot  the  biharmonic  problem  15,6)  and  the  stationary  Stokes 
problem  |10).  The  results  ol  11,21  have  also  been  used  to  analyze  a variety  ol  variation- 
ally  posed  piohlems  that  are  not  of  foim  (l.l). 

Them  ale  other  problems  of  a similai  nature,  however,  where  attempts  at  using  the 
ideas  of  11,2,51  weie  not  entiiely  successful  since  not  all  of  the  abstract  hy|v>t  hoses 
were  satisfied:  specifically  the  hiezri  condition  (1.0  or.  equivalent  ly,  the  HahuSka 
condition  (1.(0,  (s  not  satisfied  with  the  usual  choice  of  norms,  i.o.,  the  approximation 
met  tuxls  foi  these  piobloms  ate  not  stable  with  respect  to  the  usual  norms.  This  is  the 
ease,  fot  example,  in  the  analysis  in  (71  of  the  lleimann  Miyoshi  111,14,17)  mixed  method 
fot  the  bihatmonio  problem.  In  the  analysis  of  this  method  a natural  choice  foi  both 


II  •!  and  II  •'  t the  l'  ordot  Sobolev  noimt  howevei  this  method  is  not  stable  with 

V w 


respect  to  this  choice.  As  a result  ol  tills  difficulty,  the  el  1 ol  estimates  obtained  in 


n 


a 


j. 


[71  are  not  quasi -opt imal . A similar  difficulty  arises  in  the  analysis  of  the  Hermann- 
Johnson  [13,14,15]  and  Ciarlet-Raviart  [9]  mixed  methods  for  the  biharmonic  problem.  In 
later  work  of  Scholz  [211  and  Rannacher  [19]  quasi-opt imal  error  estimates  were  obtained 
for  the  mixed  methods  of  Ciarlet-Raviart  and  Hermann-Miyoshi  , although  the  systematic 
approach  of  Brezzi  and  BabuSka  was  abandoned. 

In  a forthcoming  paper  of  BabuSka,  Osborn,  and  PitkA'ranta  [31  quasi-opt  imal  error 

estimates  for  mixed  methods  for  the  biharmonic  problem  are  derived  by  an  application  of 

the  results  of  Prezzi  and  BabuSka.  In  this  work  a new  family  of  (mesh  dependent)  norms  are 

introduced  with  respect  to  which  the  above  mentioned  mixed  methods  (Ciarlet-Raviart, 

Hermann-Miyoshi , Hermann-Johnson)  are  stable.  Error  estimates  in  these  norms  then  follow 

directly  from  the  results  of  Brezzi  and  BabuSka,  once  the  approximation  properties  of  the 

subspaces  V.  and  W,  have  been  determined  in  these  new  norms.  Error  estimates  in  the 
h h 

more  standard  norms  are  then  obtained  by  using  the  usual  duality  argument. 

It  is  the  intent  of  this  paper  to  provide  an  abstract  approach  to  the  analysis  of 
mixed  methods  which  leads  to  quasi-opt imal  error  estimates,  uses  only  standard  norms, 
and  is  systematic.  We  shall  assume  that  existence  and  uniqueness  for  the  continuous 
(infinite  dimensional)  problem  has  been  established  and  develop  an  abstract  framework  under 
which  quasi-opt imal  error  estimates  can  be  derived  for  a variety  of  examples  which  do  not 
fit  within  the  convergence  theory  of  Brezzi  and  BabuSka  using  the  usual  norms. 

Section  2 contains  the  abstract  convergence  results  of  the  paper.  In  Section  3 we 
present  four  examples  previously  analyzed  in  the  literature  and  show  how  error  estimates 
can  be  derived  from  the  theorems  in  Section  2.  Three  of  these  methods  are  mixed  methods 
for  the  biharmonic  problem  and  the  fourth  is  a mixed  method  for  a second  order  problem 
analyzed  by  Raviart -Thomas  [20] . 

It  is  interesting  to  note  that  in  this  last  example  the  results  of  Brezzi  and  BabuSka 
apply  with  the  choice  of  spaces  used  by  Ravi art -Thomas , but  fail  to  yield  quasi-opt imal  error 
estimates  in  all  cases  due  to  the  way  in  which  the  variables  are  tied  together  in  the  error  est  i- 
mates.  In  our  analysis  the  error  estimates  for  the  two  variables  are  separated  and  quasi- 
opt  imal  error  estimates  are  obtained.  For  the  throe  mixed  methods  for  the  biharmonic 
problem  that  are  analyzed  in  Section  3 the  results  of  the  present  paper  and  those  obtained 

-3- 


in  13).  using  different  techniques,  ate  (he  same,  for  additional  lesuUs  on  nux.d 
»e«  Oden  ( 18  J . 

m , p 

w * 


Throughout  this  paper,  we  shall  use  the  Sobolev  spaces  w"’1'),.),  wher. 
convex  polygon  in  the  plane,  m is  a nonnegative  integer,  and  1 . j-  . 
we  have  the  sominorms  and  norms 

1 ' m.p.Q  ’ 1 [ , J lc''‘vlr  'utl  P 

and 


i a 


On  thes«  s j . 1 1 


« a I 


11 ra , p , il  -</.,!  lpVr 

i*  » * m 


dx> 


1 l' 


When  j- 


and 


we  denote  w" ' * ( ,1)  bv  h"  ( il)  and  write 


v „ , - v 

n> , m , »: 


»v«  - Hvll 

® ^ * **  TO  , >1 

We  will  further  denote  by  nJ^Oll  the  subspace  of  W1'^:)  of  functions  that  vams 

" • “ 1 and  **  "o‘«>  th<‘  suhspace  of  of  functions  that  vanish  together  with 

thair  normal  derivatives  on  r.  For  m-1  and  2 we  will  also  use  the  spaces 

-m  id 

H «*;>  - IH0 <»!)1*  (the  dual  space  of  h"(,s>)  with  the  norm  on  taken  to  be  the 

usual  dual  norm.  To  further  simplify  notation  we  often  drop  the  use  of  the  subscript  >: 
m the  norm  when  the  context  is  clear. 
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Al'stt.-nt  Ko  su  Its 


II  , 1 t!  . .tint 
V W 


hlu.v'  ! - «b«  !'ullv  II  v I! 


lot  v,  w.  and  H bo  t hi  oo  real  ttanach  *pxtvos  with  n turns 
i osp«n  t ivo l y . Wo  annum*  V » H with  .t  cont  mucus  imbedding  . l ot  .1 1 • , • ) .tin!  M* 
N’  oont  inuoutf  bilinear  forms  on  H \ H and  V > W , respectively. 

(-.1)  aUi,v)  | ^ Hall  Hull  Bvll  V u,v  « M , 

v u c v,  v v * w . 

Wo  vonsidei  t ho  following  problem,  which  wo  tofet  tv'  ,v»  problem  V 
r;iv«»t\  f V and  g » W , find  (u,v)  c V * W satisfying: 

(-•*'  .t  (u,v)  * biv , , ' - (f  ,v>  V v .•  v , 

(*  4'  Mu.o>  - Ig ,/)  V 0 < W , 

whc'i  o (•»•'  denotes  t ho  jvtiting  bt»t  woon  \ and  V 1 01  w and  w* 

"v  »h*U  ho  vnt.otosf.Hl  in  this  proMom  fov  .1  of  t.o  toi  u..,v 

who  1 o t'  is  a subcUss  of  V v W.  Wo  shall  a* sumo  that-. 

UU'  Kot  if. o'  1',  I'  has  a uni.juo  solution. 

In  tho  AHAlysis  of  I'toblnm  y wo  will  Also  vonsi.loi  tho  A.lvomt  Mom 

■'  • *'»'«’« o u vs  A banaoh  sjmoo  satisfying  w with  .1  ,-ont  lvuvou ■■ 

u-.vNM.Hng.  f.ml  iy.il  . ly^l  v x M sAtvstvvng 
1 •'  • s ' a 1 v , y ' • h ( v , V ' - 0 V v 1 V 

- («!,.•'  , v »■  . w , 

Wo  shAll  ASSIURO  t IlAt 

iM.  ' rioMom  v.'.t.v  h«*  a uvu.juo  solut  von  fot  OaoIv  .1  . • . 

rhioughout  this  papoi  wo  shall  bo  oun.otn.xl  with  tho  problem  ot  ,t|1>iox,*.it  uw  th. 
I vv  t in,.  ot  l iowAi.l  this  on.l , wo  suppose  wo  aio  givon  fmito  .1  ivvvon  ■ ions 

SI-A.O-.  Vh  \ ttv.l  «h  K.  Wo  thon  oonsi.loi  tho  following  appioxiniAt  o ot  oil  .0  ... 

v...  lofoi  to  As  oiol'lom  V 

h 

I'm. I ' V?1  ' Wjx  SAtisfying: 

'••  •'  Al«„.v'  . blv.,hi  - V f , v ' V V . vh  . 

'*  fi'  b(U  • lv),v'  ' V V*  W 

•'  h 

w«*  Will  thon  Viow  .IS  an  appi ox imat tor  tv'  u and  as  an  appioxiaat  u>n  ti'  v 

In  thin  si-i't  ton  wo  '.tin  ost  im.itos  tv'i  u u and 

h * *h  • 


_ s 


I 


We  now  stato  several  further  assumpt ions  which  we  will  require  in  the  proofs  of  our 


main  results. 


There  is  a constant  i N 0 ( u independent  of  h ) such  that 

2 

a(v,v)  ■ illvlP  V v . z,  , 


wht'tf  2 tv  V : b(v,»)  - 0 V . W ) . 

t.  n 1) 


S(h)  is  a number  satisfying 


llvll v _<  S(h)  llvll 


V ' Vh  • 


There  is  an  operator  n.  : Y * V satisfying 


b(y  - thy ,v)  = 0 


y < Y and  V yf  t W, 


where  Y - span  (( y^  ^ , , u)  , (u,i>)  is  the  solution  of  problem  P , and 

(y^,  \ ,)  is  tln>  solution  of  (2.5)  — (2.6)  corresponding  to  d . G*. 


For  the  examples  treated  in  Section  3 the  existence  and  uniqueness  of  the  approximate 
solution  (u^,i^)  carl  ^e  established  in  various  ways.  We  now  give  a proof  based  on  the 
assumpt ions  made  above . 


Theorem  1.  Assume  that  hypotheses  (H2)  , (H J)  , and  (H5)  are  valid.  Then  problem  P 

h 

has  a unique  solution. 

Proof . Since  and  are  finite  dimensional,  it  suffices  to  show  that  if 

(u  , v.  ) • V ' W satisfies 
h h h h 


a (u,  ,v)  + b(v,  ) = 0 
h h 


V ' Vh  ' 


b(u  ,>f)  = 0 
h 


V*  * , 

h 


then  uh  = = 0.  Choosing  v = u^  in  (2.9)  and  *•  = — 0 in  (2.10)  and  adding  the 

equations,  we  get  a(u.  ,u.  ) *>  0 . 

h h 

Noting  from  (2.10)  that  u Z and  using  (H3)  we  have  llu  II  = 0.  Hence  u = 0 . 

h h h H h 

Setting  u.  = 0 in  (2.9)  we  obtain 
h 


b(v,ii-  ) - 0 V 

h 

<d,<  ) 

**h*e  TSTT, 

diC  G* 


V 1 Vh  • 
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Bv  tH2),  for  each  d « d',  there  exists  yd  . v such  that  for  all  *-  t w 


(d.c)  - b(y 


Thus 


(d'V  - b(yd,V 


' b(Vd'*h 


(applying  HS) 
(using  (2.  111). 


Equation  ( - . 12)  then  implies  ■*>  - 0 . 

h 


Oui  main  result  in  this  section  are  Theorems  2 and  3 which  present  abstract  estimates 

h - 


for  the  errors  u-u,  and  < -o 
h 


Theorem  2.  Suppose  hypotheses  (H1)-(H5)  are  valid  and  that  (u,H)  and  (u  . tf.  ) are  the 

h h 


respective  solutions  of  problems  P and  P^  Then  (with  it  defined  by  (HS)) 


<2.  13) 
and 
<2.  14) 


l,U~Uh"ll  — a + (Hall  + a)  Hu-n^ul^l  for  a11  *'  ’ wt] 


S(h) 


«“-uh»v  i «U-V"v  + — a — III  blls  (h)ll  ii>— V II  w + Hall  llu-tt^ull^)  for  all  . 

I!  in  addition 


(Ht>) 
who  1 1 


{v  . V:  b(v,v')  - 0 V v'  r W}, 


t hon 
UM5) 
and 
C.  lb) 


"U-Uh"H  - 


1 + 


Hall 


llu-n.  ull 
)t  H 


"-Vv  i IIU-nhU(V  + “I—  llu-”hUllH 


PrvXJl  . 

(2.171 


s inq  (2..))  wo  see  that 

a ( n^u , V ) + b(v , i)')  - a (u , v)  + b(v,t(i)  + a(tt  u-u,v) 


(f,v)  + a(*  u-u,v)  V v , V. 

h ti 


and  ft.m  (.’  4)  and  (HM  we  see  that 

<•’.181  b(s  u,c)  »(q,v")  V 

h 


f W. 


subtracting  (2.7)  from  (2.17)  we  find 

( - • 1 a(t  u-u  ,v)  ♦ b ( v . ti—  !;■  ) » a( 
h h h 

and  ••‘iiht  t icting  (2.9)  from  (2.18)  we  obtain 


V v , V 


h ' 


1 


....  


-7- 


(2.20) 


b(  ”,  u-u,  , »■)  - 0 V 

h h 


“h  ‘ 


Choc-.ing  v - u-u  in  (2.14)  we  lidvi 

h 


a (”,  u-u.  , ”,  u-u  ) * b(”,  u-u.  , 4— »,  ) - a ( < u-u , * u-u.  1 . 

!i  h h h h h h h h 

Apj -lying  (2.20)  we  get 

>0.21)  .11  •>  u-u  , ” u-u,  ) - a(”  u-ut«  u-u.  ) * b(u, -j.u.i-v"  ) foi  all  0 w 

n h h h h h !>  h h 1 


Osinv)  (2.1),  (2.2),  (H3)  , (!14)  , ami  noting  from  (2.20)  that  ••[u-ul  , 2,,  wcthencbtai 

ili  i^u-u^r  ~ ^ Hall  II  i^u-ull  (|ll  n^u-u^ll  H ♦ KbBsthlllu^-iij  ull  U 0-vllw 

and  hence 


(2.22)  Ili.u-u  II  - [II  all  II  u-”  ull,  ♦ D bU  S (h)  H *->'11  ) for  all  e . W 

h h H — a h H W h 


1 , 


Thus 

II  u-u.  I!  : u-e  ull  * H”.  u-u  II  - —III  Wl  S (li)  II  v — 'll  * (II  all  ♦ a)  llu-”.  ull  1 

h H h H h h H — a W h H 

for  all  c t w,  This  proves  (2.1)). 

In  order  to  prove  (2.14)  we  first  note  that 


lu-u  I - Hu-”  ull  ♦ lh  u-u  II  v |lu-”, ull,  ♦ S (h)  ll  - u-u  II 
h V — h V li  h V — h V h h 11 

(2.14)  now  follows  from  (2.22). 


To  prove  (2.1b)  we  observe  that  (2.20)  together  with  2,  • 2 implies  that 

li 


(2.21) 


b(  ".u-u,  ,»•)  - 0 
h h 


Hence  (2.21)  simplifies  t o 

(2.24)  a ( ",  u-u  , ” U-u  ) ■ a ( n.  u-u , n.u-u.  ) . 

h h li  h h h h 

Applying  (2.1)  and  (H>)  to  (2.24)  yields 

(2.2b)  Kw.u-u.ll  v R n u-uH 

h h H — a h 11 

(2.1b)  follows  by  the  triangle  inequality. 

Ti'  establish  (2.1b)  we  wiite 

II  u-u  H • l!u-n  ull,  » II”, u-u,  II,. 

h \ - ll  V h h V 

Hu-n  ull  * S (h)  ||  - u-u.  II  (by  114' 
ll  \ h h H 

; " u — ull  * S (h)  |l u-”  ull 


(using  (2 . 2b) ' . 


Gorol  lai  y 


Inequality  (-.15)  holds  without  assumption  (H4)  . 


Theort-ff \ . 

a'  Suppose  hypotheses  (Hi),  (H2)  , (H3)  , and  (H5)  are  valid  and  that  (u,„)  and  (u^  , . }) 
are  the  respect ive  solutions  of  problems  P and  P . Then  (with  (y  , X ) and  as 


d d 


defined  in  (H2)  and  (H5) , respectively)  , 

U.2€0 


'♦■♦h'c  ’ ^P(b<Yd-"hV<-'>  * alVU'Vd'V  4 btu-uh'\ 

at  G 


-n)  VI!  d! 


for  ail  * , n < w 


b)  If  in  addition  (Hi)  holds  (2.  *•  Z)  , then 

h 


(-'.-'7) 


b(u-u.  -n)  b(u-*  u , X ,-n) 

h d h d 


n * W 


c)  If  we  further  have  that 

(H7)  There  is  an  operator  E : A * wv  satisfying  b(v,}\  \-\)  *=  0 for  all  v V and 

h h h h 

^ -V  where  A * span  ( { .,,;),  (u,v)  is  the  solution  of  problem  \ , 

and  (Vj#\j)  is  the  solution  of  (2. 5)-(2.t>)  corresponding  to  d G'  , then 


(-'.-'8) 
and 
(2.2*) 

Proof.  From  (-.6)  we  have 
(2.30) 


b(yd"Vd'v'!h<‘)  ‘ (d'v_;h*’) 


b(u-VVW  - ‘«'Vrhxd’  • 


II.-*  II  . - sup  (d,  c-v  ) /II  dll  . , - sup  b(y.,;-v  ) 
dt  G 1 h ° dcG*  d h 

Subtraction  of  (2.7)  from  (2.3)  and  (2.8)  from  (2.4)  yields 

(2.  ill  a(u-u.  ,v)  * b (v , ii-i,  ) 0 V v t V 

n h 

and 

(2.32)  b(u-u  , n)  ”0 

II 


II  dll 


h 


combining  (2.S),  (HS) 

. (2. 

31),  i 

ind  ( 2 . 

32)  we 

obtain 

b(V**V  ■ 

blyd 

_V’d‘ 

■*-v 

4 b(v. 

m 

b<yd 

■Vd- 

•4-v* ) + 

a(uh‘U 

• Va* 

• 

b(yd 

■Vd’ 

„ -v' ) ♦ 

a (u  -u 
n 

■vvyd' 

♦ a(u.  -u, 
h 

YV 

* 

b(yd 

-Vd- 

. -V ) ♦ 

a (u,  -u 
h 

■vw 

♦ b(u-u,  , 
h 

V 

b<yd 

hYd  ’ 

V’"V)  4 

a ( u -u , 
h 

’W'V 

♦ b(u-u.  , 
h 

\ - 
d 

Subst itut ion  of  this  identity  in  (2.30)  yields  (2.26). 


If  Z «-  Z then 
n 

b(n  u-u  ,*p)  - 0 V «£  . W (see  (2.23)  above) 
h h 

and  so  (2.?7)  follows  immediately. 

Now,  if  in  addition  (H7)  holds,  then 

b(Yd_7Thyd#  ^”‘‘h^^  = bfy^-^Y)  = (d»^~E.0)  (by  (2.6)). 

and 

b(u-u  , A -E.A  ,)  = b(u,  \-)'  A.)  = (g,^-E.X.)  . 
h d h d dhd  ^ d h d 

Thus  (2.28)  and  (2.29)  are  established. 

Remark.  Note  that  inequality  (2.15)  in  Theorem  2 and  all  the  results  of  Theorem  3 hold 
without  assumption  (2.2).  This  observation  is  used  in  subsection  3c. 

We  end  this  section  with  several  remarks  on  the  hypotheses  (H3)-(H7).  We  assume  here 
that  V and  W are  Hilbert  spaces. 

1)  It  is  clear  that  if 

2 

(2.33)  a(v,v)  > yn  II  v||  for  all  v c , 

— 0 V h 

then  hypotheses  (1.3)  in  Brezzi’s  theorem  is  valid.  In  the  applications  we  consider  in 

section  3,  (2.33)  is  not  true  (with  Yq  independent  of  h ) but  is  valid  when  llvll^ 

is  replaced  by  llvll  . This  accounts  for  (H3)  (and  (H4)). 

H 

2)  In  hypotheses  (H5)-(H7)  it  appears  that  we  are  not  making  use  of  conditions 

similar  to  (1.4) . In  fact,  in  applications  the  operator  described  in  (H5)  is  often 

constructed  in  order  to  verify  (1.4).  A more  precise  relationship  is  given  below  in 
Propositions  1 and  2.  For  further  ideas  in  this  direction,  consult  the  work  of  Fortin 
[11]  . 

Proposition  1 . suppose 

(2.34)  sup 
V'Vh 

where  k > 0 . Then  there  is  an  operator  it  s V ► V that  satisfies 
0 h h 

b(v-u  v,^)  =0  V v t V and  V < W 

h n 

and 


|b(v,l) 


i koIMIw 


r 

k’ 


) 
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V 


V 


IlK.Vll  • ~ II V# 
h V k0  V 


V . 


V,oof  v'onsol.'i  h(v,  »•}  on  ?..*  «W.  .where  2.‘  - |v  , V.,  v ia  V-or  t hoqona  1 to  1 


h h 

Wo  i mintM  i at e ly  see  that 

sup  N kJMIM  V **  f W4 


ami 


i N VII..  0 W 

v,  Zh  V 


»up  |b(v,v")  | ' 0 V 0 ^ V < Z 


w, 


h 


It  thus  follows  from  |1,2|  that  for  each  v . V t hei  e is  a unique  *hv  . 7 satisfying 

bin  v.y'l  - b(v,v’)  V v'  » W 

h h 

Kill  t hoimoi  e , 

ti  i,  Ill’ll  „ n 
lln.vll  <•  — r — llvll.  . 

h V ~ V 

This  proves  Pioposition  1. 

Wt>  Also  note  t h«*t  it  follows  from  (1,21  that  for  ear  h ^ » w t her*'  exists  a 

unique  ' , v*  W.  satisfyinq 
h h 

Mv,'  v'>  - I'Vv.vM  V v , 2l  . 
h n 


Pm  t hot  mot  o , 


« ! . -•«  111,11  IU-II 

li  W W 


Pi  op'-i  i t ion  . Suppose 


1*’ »'l  1.  ..  K||^n 

1 (Fvfl  i 

v.  V V 


V 0 * W , 


whew  k 0 , ami  sup\xvu'  there  is  an  operator  n tV  * satisfying 


b(v-\v,v’)  - 0 v ¥*  « w. 
h n 


In  vK  • Oil  vll  V v » V . 

h V V 


Then  ( .* . 141  hoi  its  . 


II 


\ . Applications 

In  this  section  wo  apply  the  results  of  Section  2 to  several  examplo  . 
a ) v j at lot-Kaviart  method 

Considet  the  bi harmonic  problem 


b-  2 


0 on  r - DU, 


where  U is  a convex  polygon  ami  g is  a given  function.  If  g i I!  ( ) then  t hi  i • i 

a unique  solution  i*'  » of  (1.1).  In  addition  the  following  rcqulatity  n nil  \ 

- 1 3 2 

known  for  this  problem:  If  g « H (ft)  , then  . H (ft)  n H^(ft)  and  thoie  i a 

constant  C such  that 

(1.2)  II  ij'll  ^ jc  ell  gll  ^ V g e H (ii)  . 

Using  the  well-known  correspondence  between  the  biharmonic  problem  and  the  stoke-,  pioblem, 
this  regularity  result  can  be  deduced  from  the  regularity  result  for  the  Stokes  pioblom 
proved  in  (161  (cf.  also  (121). 

We  now  seek  an  approximation  to  the  solution  g-  of  (i.l)  by  a mixed  method,  i.e.t 
we  introduce  an  auxiliary  variable  (u  -Atf*  for  the  method  of  this  subsection!  wiite  (\  l) 
as  a lower  order  system,  cast  this  system  in  variational  form,  and  then  con  \dci  t h«  kit 
Galerkin  method  corresponding  to  this  variational  foimulat  ion.  In  part  ieula»  , the  mixed 

method  we  study  will  be  based  on  the  following  variational  formulation  ol  (t.11,  « it  :-t 

considered  by  Piarlet  and  Raviart  (91: 

Given  g « H 1 (ft)  , find  (u,i|>)  « ll'  (ft)  x H*(ft)  satisfying 


/ UV  dx  - / Vv  • V \J<  dx 

ft  ft 


V V t H (ft) 


* / Vu  • W dx  » -/  gv'  dx  V v*  ■ H 1 ( >:>  . 


Using  the  regularity  result  (t.2)  it  is  not  difficult  to  show  Theorem  1 ol  (“1' 

that  if  tj/  is  the  solution  to  (1.1)  and  u - — Ai|» , then  (u,tj)  is  a -.elution  of  (t.l), 
ami  if  (u,iji)  is  a solution  of  (t.i),  then  ^ is  a solution  of  (l.l)  and  u b 
It  is  clear  that  (1.  0 Is  an  example  of  problem  p of  Sect  ion  .*  with  V II*  (s')  , 


W - H^(ft),  H I. , ( il)  , a(u,v)  - f uv  dx , and  b(u,\j)  -j  \\i  dx, 


(ami  with  g t »-pl  a*  ■ .t 


by  -g)  . Here  the  subclass  P of  data  for  which  (Hi)  is  satisfied  \ given  b' 


-l  1- 


D * O * W*.  Since  the  form  a is  symmetric,  the  adjoint  problem  (2.5),  (2.6),  with 
G » w » H^(i.)  , is  the  sane  as  problem  P and  thus  is  uniquely  solvable  for  all  d » W . 
Hence  (H2)  is  satisfied.  Usinq  (3.2)  we  also  have 
(3.4)  !lyd!l  «■  I!  \dll3  £ Clldll  . 

Next  we  discuss  the  finite  dimensional  subspaces  used  in  the  approximation  scheme. 

For  0 < h • 1 , let  t ^ be  a tr ianqulat ion  of  ft  with  triangles  T of  diameter  less 

than  or  equal  to  h . We  assume  the  family  {x^  } satisfies  the  minimal  angle  condition, 

i.e.,  there  is  a constant  o % 0 such  that 
h 

T 

max  — <o  V h , 

P ~ 

Tt  T.  T 
h 

where  h^  is  the  diameter  of  T and  p ^ is  the  diameter  of  the  largest  circle  contained 
in  T,  and  is  quasi-uniform,  i.e.,  there  is  a constant  x > 0 such  that 


— — — < x V h . 

min  h — 

_ T 
T 

For  k 1 a fixed  integer  we  define 

(3.5)  Sh  = (v  i C°  (jj)  :v  |T  c Pk  V T e Th  } 

where  P^  is  the  space  of  polynomials  of  degree  k or  less  in  the  variables  x^  and  x., . 
We  then  consider  the  approximate  problem  P^  with  and  n H^(S1)  . Note 

that  this  scheme  yields  direct  approximations  to  iJj  and  u = -A1J1  (the  stream  function  and 
vorticity  in  hydrodynamical  problems) . 

To  apply  our  theorems  we  must  check  that  hypotheses  (H3)-(H5)  are  valid.  (H3)  is 
clearly  valid  with  a=l  and  since  our  family  of  triangulations  is  quasi-uniform,  (H4)  is 
satisfied  with  S(h)  = C/h  for  some  constant  C . It  remains  to  check  (H5) . For 


v f H (S!)  define  x.v  by: 

h 


V ' V 


J Vfir^v)  d x = / Vv*w>  dx  V t V 

n si 


j xv  dx  = / v dx  , i.e., 

S2  " Si 


Then  (MS)  is  sat  i sf  i*-d,  and  in 


« v l*e  t he  Neumann  projection  of  v into  V . 
h h 


•dil it  ion  »t andard  approx imabil lty  result*  imply  that  it  v • H '(>•),  i 


(i.t>)  II  v - w vll  • Ch  II  vll  , 1-0,1  and  1 » i min  ( k*  1 , r-2)  . 

h ) — ( 


We  are  now  ready  to  apply  Theorems  2 and  3.  Suppose  * H (>•'),  i 


k 2 . Then,  using  (2.13),  (1.6),  and  standard  approx  imahi  1 1 1 y i exult  a, 


( 3.  /) 


-1 


lu-u  1 C (h  inf  M , ♦ Hu-" , ull  ) 
h 0 - . 1 h 0 

h 


• C (h” 1 h*  " lll  it’ll  » hM**llull  ’ 

S S-i 


• Ch*  “II  ^11  (since  u - 

s 

where  x * m m (r , k«  l ) . 

From  (2.14)  we  find  in  a similar  fashion  that 

( 1.8) 

whet e s min(r,kt l) . 

Finally,  from  (2.1),  (2.2),  (2.26),  (1.4),  (l.7\  and  (1.8),  we  have 


JIu-uJI.  Ch  ' 'll  ^||  , 

hi—  s 


■ ‘'““Ci,  J*vd‘*hvdli  inl  lu'^"1 

d»H  (O  *tW. 

h 


* '“•"h'o'Vd'Vd'o  ' ""-Vi  ,,,f  11  Y'V/,M"-i 


ni  w 


1 


I-  am.  Illy  ,11,  h ' II  <41  t h ‘11*11  h II y II 

d.H-V>  d 1 “ M d 


s-  « 


. h iuii  ir  ii  \ ii  , i/imii  , 

a .1  ' -1 


C h8‘l|li|.|l 


whet  e mi  n ( t , k ♦ 1 ) • 


Since  ( t .7)  - ( t . 4)  are  valid  only  for  k 2 , the  met  hods  of  this  papet  do  not  yield 
error  estimates  for  the  case  k • l in  this  example.  For  this  case  the  i oadet  is  referred 
to  Scholr  |22].  The  estimates  (1.7) -(1.4)  improve  on  those  in  Ciat let  and  Kaviart  (4| . 
Scholr  1211  obtained  (1.7)  under  the  assumption  that  1*  is  smooth.  (1.7)  and  (1.4) 
were  obtained  by  BabtiSXa,  Osborn,  and  Pitkaranta  (l). 
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1 


) 


J 


w«*  t *«ia  i k that  Th«*oi  t'«n  ' could  a l no  be  used  to  obtain  .in  on  or  vat  im.it  «•  for 

4 

II  ||  (by  . boosing  Cl  - l, (ill)  when  ii  • " <>:)•  M' 'wove i in  ordoi  to  got  gu.i; i -opt  im.i  1 
h 0 

4 


1 1 


Sults  wo  would  require  t ho  regularity  result  that  A * \ %(iJ)  implies 


A 


H ( ) , 


• huh  is  not  valid  on  a convex  poison. 


b)  Hermann  Mi^oshi  met  hod 

Wo  consider  in  this  subsoot  ion  another  mixed  method  for  the  approx  imat  «•  -elution  *m 
(i.l).  in  this  method  the  auxiliary  variable  is  the  vector  of  second  partial  derivative: 


l ot 


V - tv  - (V . p . 1 • i,i  ' - ! Vj,  - v2l, 

(with  thf  usual  product  mum),  .tn.l 
1 , 

W - It  lill  . 

0 

Thou  tho  ml  x ml  mo  tin'll  wo  study  will  l'o  based  on  t ho  following  variational  lonnul.it  ion 
ot  l t.  II  . 

-1 


t <■  li'l 


ill  von  i)  . H l(ii),  find  (u.i|l  . V v w satisfying 

| V j u v dx  t ' / T— 'lx  » 0 

j i/i-i  L " ” i.V-i  >:  ''xi  ?xi 

\ 


V V r V 


V f ~~r^~  dx  - - f UV  dx  V 
,'x  ,'x 

1-1  i.'  t i >: 


Its  i mi  tho  rouulat  ity  result  O.J)  it  is  not  difficult  to  show  that  it  <•  is  a solution 

of  (t.ll  and  yi  - (u  i«  defined  by  u » — ' , then  (yi.il)  is  a solution  of 

( l . I o > . and  it  t\i.vl  is  a solution  of  (.1.10).  then  ij  is  a solution  of  It. 11  and 

U * " , 

ll  dxpx 

Wo  oasily  observe  that  (t.lt'l  is  an  example  of  problem  P with  V and  w as  above, 
It  - iv  - (v,.),  1 • i,t  - - :V  -v  , v < I.  ,(lD) 

V 11  — l v . 1 t .1 

Iwith  tho  usual  product  norm), 


,i(ii,v)  V f u v dx  , and  b(u,i)l  - ^ I 

4 ' i/i-l  i)  11  li  ' i.'i  i »: 


■hi 

3X  ,'x 
) i 


dx  . 
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As  in  Subject  ion  .'a  the  subclass  D of  data  fox  which  (Hi'  is  satisfied  » .i\,  . * 

p • 0 \ v*'  , and  since  a is  again  symmetric  (H.M  is  satisfied  with  c.  w h‘»  ' 

letting  be  as  defined  in  O.M  , we  then  considei  the  apptoximat.  ptobltw  », 

h 

with 


and 


Vh  ’ ‘t  * (vij>  : V12  - V21  • Vi.1  * V 


wv  * S.  a H_(Q)  . 
h h 0 


With  this  choice  for  the  forms  a and  b and  the  spaces  V.  and  W,  , pi  obi  cm  v 

h h h 

now  describes  the  Hermann-Miyoshi  method  jO,  14,  1 7 1 for  the  approximations!  t'.  1 ■*.:  nr  •••■.» 
ptoMem.  Note  that  with  this  method  we  obtain  direct  approx imat  ions  to  w and 

x .*x 

i 1 

(the  displacement  and  the  moments  in  elasticity  problems'. 

As  in  Subsection  'a  we  have  hypothesis  (HO  satisfied  with  . % - 1 and  iH4'  -it  i * i«  «i 

with  S ( h)  * C/h  for  some  constant  C . (H5»)  for  this  example  is  contained  in  lemma  ' 

in  (7).  Mot  cover,  by  a minor  modification  of  the  proof  of  lemma  in  171  we  obt  un  tv, 

existence  of  * : V k V,  satisfying: 

h h 


b ( v 


v , C ) - 0 


r-2  4 

and  for  v , V > [H  (».'))  , r - the  estimate 


('.11'  II v - «*  vll  • chV  Slyll 

' n i \ i 1 


1 ~ » 1 And  1 ( • min  (k«  1 ,1  -2)  . 


Wo  oan  now  Apply  Theorems  2 And  ^ in  the  same  way  as  in  ibsoot  ion  'a  .\vP-iiun.i 

thoso  thooroms  with  (l.ll)  And  stAndard  Approx i mAbi 1 1 ty  results,  wo  obtain  f.  i n i 

with  t ^ 3 And  k ^ 2i 

0.12)  II u-u  |l  . o h"  ' II 

' vn  U ***  •; 

F - V 


II V”  vh11 1 - c h 11  ill , 


II  *-<^H  J » i'  h R 'll 


o.i  n 

and 
0.14' 

where  s - nin(r,krl)  . Estimates  (i. 121-0. U)  improve  thoso  initio.  ■ Knvi.ut  '1 
KAniiAoher  11^1  rooently  proved  these  estimates  foi  k - 2.  IiAhu-'fk.i . Nm  n .v 
1’  i t k a r an  t a (<1  proved  (t.12)  and  (1.14'. 
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c)  Hermann -John -on  method 


Wi  . otisidei  here  a further  mixed  method  t ox  t he  solution  of  (3.1)  in  which  the 

auxiliaxy  vai  i the  v«  tot  v’t  t * ond  partial  derivatives  of  v # as  in  section  3b. 

Given  a triangle  r k and  a function  v - (v  . ) with  v t H*  (T)  . 1 * i.i  *•  2, 

ii  ij  - J - ' 


and  v 


ami 


i; 


:\ 


WV  drt  in*’ 


M (v)  ) v v \ 

i.j-i  11  > 1 


M (V) 

vt  ' 


y 

i,‘)-l 


v . . v. t . 

ii  \ t 


where  v (\^  v,)  is  the  unit  outward  normal  and  t » { i ^ , t , ) - (v>#-v^)  is  the  unit 
tangent  along  »T.  Let 


V V (h) 


tv  (v  ) : V . . 1.  (i:)  , v 

i i ij  2 1 


21  ' 


and 


v H (T)  v T * : , . and  M (v)  is  continuous 

i ’ T h v v 

at  the  int  ere  lenient  boundaries)  with 

ii  yC  V T ii  v.  id 

' V , ^ 1 ij  ,T 

i » i'  l t*  t J 

h 

w W / 1 u.)  , where  \ is  some  number  larger  than  2 . 

The  mixed  method  w«-  study  in  this  subsection  will  be  based  on  the  following  varia- 
t lonal  formulat  ion  of  (3.1). 

Given  g H *(G),  find  (u#i^)  r V * w satisfying 

.'v . . 


( 3.  IS) 


V 


[ u v ,dx  ♦ ^ { V 

, *.  . ; ij  ij  • •.,„?*  .>x. 

1.  ) l W T.  ,h  j-1  T t 1 


Ml 


r1*  dx  - f M (v)  ^ ds  ) ■ 0 v y « V 


3T 


3 \ 


L 


J I J / 'l*  - / M lul^d.}--/  qv'lx 

i.J-i  t ,xt  ’Xi  JT  v'  X •”  a 


The  correspondence  between  (3.1)  and  ( ' . 1 *^ ) is  the  same  as  the  correspondence  between  (.3.1) 

and  (3.10),  i.e.,  if  ; is  the  solution  of  (3.1),  then  ((  1 , tf)  is  a solution  of 

v3x  3x 

‘ ) 

(3.  is),  and  it  ((u..)»tf)  is  a solution  of  (3.ir>),  then  i-  is  the  solution  of  (3.1)  and 


JL 


t j 


‘X  >x 

i ) 
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One  easily  sees  that  (3. IS)  is  an  example  of  problem  r with  V and  W as  above, 
M as  in  Section  3b, 


ami 


a(^.v)  - 
b(^,\j)  « 


V 


1 u ■ jV,  , d* 
, ij  ij 


».i-i  a 


v 1 . I / w1  S'  ***  • / Tf  d"  } • 

T,  ih  i, j«l  T v 3T  V' 


3li 


As  in  the  previous  subsect  ions,  a subclass  P of  data  for  which  (HI)  is  satisfied  is 
ijiven  bn  P * 0 * H (il)  , and  since  a is  aqain  symmetric,  (H2)  is  satisfied  for  G « h'l.) 
1,h*  net e that  in  this  example  the  space  V ■ V(h)  depends  on  h.  Por  each  h the  form 
b(yi,;)  is  bounded  on  V(h)  » W (where  W * W1,P,  p > 2)  with  a bound  b that  depends 
on  h • In  t he  error  estimates  in  this  subsection  we  do  not  require  that  this  bound  be 
independent  of  h . Cf . the  Remark  following  Theorem  3. 


Letting  S be  as  defined  by  (3.5),  we  then  consider  the  approximate  problem  p 


with 


and 


‘t  ' v,vi3lT  • Vi  v TtV 


\ ’ *h  11  V0>  • 

Witt,  this  choice  for  the  forms  a and  b and  the  spaces  V.  and  w.  , we  have  the 

h h 

method  of  Herrmann- Johnson  (13,  14,  15]. 

As  in  the  previous  subsections,  hypothesis  (H3)  is  satisfied  with  .i  = 1 . We  now 
consider  (H5) . 

is  defined  by  the 


For  v v we  define  n.v  » V.  as  in  (7,  Section  41,  i.e., 

h h 


oondit ions 

( 


\ le' 


T. 


MJv  - nhv)q  ds 


V q , P 


k-1 


and  for  all  sides  T’  of 


and 


I |v  - (n  v)  |q  lx 
it  O h-v  ij 

and  v T . i.  . 


V q . P 


k-2 


By  l.emma  1 in  (7),  n^v  is  uniquely  determined  by  (3.1b).  Since  we  can  write 


t 
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' , „ i i 'x  .ix 

r.  i t , i- 1 r i ) 


nr 


M.l*i)  oaily  implies  (H5)  . Wo  not  o that  by  l.omna  4 of  | M wo  also  havo  tot  all 
I - J 4 

V V 1H  1 :>  1 , i ' » , that 


O.  171 


V " \"'  ' cht|i\'11  ; • 1 - f ■ mm<k,i  -2)  . 


Wt  next  observe  that  by  1 orwta  *'  of  17}, 


so  that  wo  at  o in  t ho  special  case-.  ot 


■iooi  orn.-.  .’  ami  In  ['articular  by  tin-  Corollary  t o Theorem  , (H4)  need  not  bo  at  i-. tied 

in  or  Jar  to  apply  (.’.15),  so  wo  shall  not  require  It  ) to  bo  quasi-uni  form.  Since  wo 

wi  h to  apply  Theorem  cat  t c,  wo  now  show  that  hypothesis  ( H 7 ) i ■,  satisfied.  As  m 

t!a-  proof  ot  l otnna  S ln  (7|  , foi  v , V,  ami  n , W - w',l'(i:)  wo  can  writ,- 

v h 0 


Hy.i.l  - Y v . i ax 

T.  i 1.1-1  T >X,  >X1 


1 1.  IH > 


A(T',y)i.  .Is  » ' b(a,y)iiia)  , 


1" . i 


h T' 


h 


u-u‘!t’  l*  ls  th*  set  of  all  sides  of  the  t ri-irujulat  ion  \ ,»l  vs  the  ;«*t  ot  all  vert  we- 
11  h h 

‘lt  ' .*  *M‘*  AIT1  ,v>  \ •-  a polynomial  ot  degree  less  t ban  oi  evpial  tv'  V m the  vaiiahlt 


■ t W we  now  choose  V V so  that 

h h 


«.1‘>> 


' hl.)  vi^x  c'  V n t x and  V T « : . 


( 1.20) 

a nd 


f ( u • J‘  i. ) q iis 
• . n 

T * 


iV-“> 


'l  « V ^ , a n d v T 1 


! !■«'  unique  olvahility  v't  this  system  is  easily  checked.  Note 


that  by  the  SoN'lev  imbed 


» w implies  ; , v*  (*:)  . Since  foi 


V.  we  have 
n 


vtuui  thev'iem, 

* Vtj  | 

lx  1W-1  ,'M,i  A ( T * , v 1 , p It  follows  ttv>m  I*.l8'  t hat  v ,,  as  del  mod  b> 

l I T !' 


( 3 . 19) -(3.21) , satisfies  (H7) . Furthermore,  by  a standard  application  of  tin  Bramblc- 
Hilbert  lemma  (4),  we  obtain  for  all  u > W n H*  (!:)  , 


Hu-.'.,  nil  . - Ch*"jllull„ 

n j — £ 


j = 0,1  and  1 < £ m in  (r  , k+ 1)  . 


We  are  now  ready  to  apply  Theorems  2 and  3.  Suppose  that  k > 1 and  ^ r.  Hr(ft), 
r 3.  From  (2.15)  and  (3.17)  we  obtain 

(3-23)  lrSh"oicllrv,o 

< C h\  IJj 


< C h II  ij/ll 


where  <5  = min(k,r-2). 


To  obtain  estimates  for  we  shall  apply  Theorem  3 in  several  different  ways. 

Choosing  G - H^(fl)  , sfi  = , and  n = (where  7^>  denotes  the  standard  Lagrange 

interpolant  of  ^ in  S^) , we  get  from  Theorem  3 (a) -(b)  that 


h"i  = ^P^'Xd  - VCd'  * - w 


d<H  (11) 


+ a(,eh  - ^'"h^d  - + b(»d  - V'xd  - Ihxd}/lldll-i  • 

To  estimate  the  terms  in  the  above  expression,  we  introduce  the  affine  transformation 


x = F(x)  = Bx  + b 


mapping  the  reference  triangle  T with  vertices  (0,0),  (1,0)  and  (0,1)  onto  T,  and  set 
v(x)  = B*1  ^ o F (x) (B-1) T , 


where  v = 


v, , v 

11  12 

v v 

21  22 


, and 


v?(x)  = ^ o F (x)  . 

Using  the  standard  change  of  variables  argument,  we  have  that  if  v £ [H£(T)]4  and 
t,  v 

* c H (T) , where  1 £ £ and  2 < t , then 


2 3v 

/ I 5^  t*-  dx  = |det  b|  / y Ilii  *5- 
T i.J-1  3xj  3xi  i i.5-1  ^ 3xi 


where 


3T 


3t 


t ;v.  ;(iis  j; 

t;  ** 


* t>  i • v r*.  ds 

T.  3,113  " dCn°te  the  tanger,t  and  1 

1 T.w  ..  , ^ ^ nonnaJ  tQ 
Sin,.  „m  , , ‘i  ' le*W  of  T.  . ' resPective 

- Ch'  '***  »!  < Ch2  (lB-i„ 

follows  that  - °/h.  and  lx.| 

' V - c>>  (c  f.  f8JJ 

" < lt:  easi 
3v 


0.24) 


if  l 

T i.j-1 


JJ  .V 


,xi  ?x 
J 1 


dx 


/ M (v)  ilf.  . i 
3T  »i  v ds) 


3r 


>v 


- Ch2f  / / x 

T i.j-1  3x.j  3x 


ii  a 


*7  :d*  + / ( f v2  ,i/~Vv;  . . 

3T  i,j  = J ij  ds} 


1Ch2"*ll1,illf»2/r 

N°W  fro,n  -1  of  ,7, ' 


0.25) 


and 


W-S-o  « v,  - . . 

1 


* V 


Slng  the  — P-Perties  of  thv.'V 

(3.27)  a . ht  lntorpolant  (c  f 

V • 0 if  y:  , p f8))  We  3iso  gpt 

ind  ky  ’ 

3-28)  ...  . . 

7f,‘  f - cerjiisf'/i 

’sing  0.24)-(j  ,o,  2>t  ‘ 

“ W°  ‘’dsiJV  obtain 


i / r _2_  , a 

- 't-vt’ij  sr-„.  v,d> 


■t  Mv> (v  " -h?’  57  <*-  - :h.-)ds 


- ch2  i 

-,Jlnf  II V - nil 

P ' rk_)  IT)  1 1,T  lnf  llv-  - „ 

1 * ' Pk(f)  q,2.f 


i Ch‘ 


e^inf  llv  - ,,i| 

r ' P .IT)  V «.T  . Il<-  - q,| 

fc-1  q e P (T,  t,T 


< Ch“ 

!^£,T  l^lt.T 

for  1 < 1 ' 

k and 

- 1 1 JL  *+l.  Changing  back  to  the  original  variables  we  further 

obtain 

Wi.f 

- Cht’3|tll.T 

and 

1'*  t , T 

- Cht"1|*lt,T  . 

Hence  if  v t 

tH  * (fi)  ] 

4 t 

and  * t H (fi)  for  1 < j,  v k and 

2 < t <_  k+ 1 , then 

'b(^  - 

\x-  * - V” 1 

(3.29) 

1 I 1 

T 

2 

^ i<%1  it  {x  - Wa  w.1*  * V)dx 

- J V(*  ' V>  57  - V)ds| 

- 1 

T 

l*ltfT  klt>T 

1 cht+e"2  llvlltHllt  . 

Choosing 

• * = \li,  4 = 1 and  t = min(r,k+l) 

s in  (3.29)  we  get 

lb(*d  ' 

\Xd'  * - V’l  1 ChS‘1|l^dlllll*lls  • 

If  k ^ 2 

we  choose  v = yj  , *:  = \ t « s-2,  and  t = 

3 in  (3.29)  to  obtain 

b(u  - 

nh“'  *d  " rhV  1 

< 

chS_1 >¥s-2  11 V 3 

< 

ChS_1  II  *•«  IIX  II  . 

5 d 3 

If  k = I we 

choose 

4=1  and  t = 2 to  obtain 

tt  u, 

, X • 

- I X 

h^ 

d 

h d 

Ch 

Mutt, 

II  X .11 , 

^ 1 

d 2 

Ch 

H U II 3 

*»d*3 

-23- 


Finally  if  k > 2,  from  (2.1)  , (2.15)  , and  (3.17)  with  1 = 1 and  l = s-2 

la<*h  - >4'  Vd  - i cll>4h  ' «"oIUh^  ' 

-cllVd  - %n0  IUhXd  - ^ll0 

< chs-2ll!ells.2  hllfcll, 

i a*8"1  ll*"s  11^!  • 

If  k = 1 we  choose  1=1  in  (3.17)  to  get 

|a()4h  ■ «»  \ ^d  ’ ja>  -ch  11  >4"  i h "y^i 

< Ch2  11*11,  II  y ,11 , . 

— 3 %u  1 

Applying  the  regularity  result 

(3.30)  II y ,11 . + II  A ,11  Clldll  _ 

%d  1 d 3 — -1 

and  collecting  terms,  we  get 

(3*31)  . — V^H  y i_  Ch  II  for  k >_  2 where  s = min  (r , k+1) 


(3.32)  II  ^“^^11^  ^ Ch  II  i|/ll  for  k = 1 . 

We  now  derive  estimates  in  1^(12).  First  consider  the  case  when  k = 1 and  * c H4(52) 
Using  Theorem  3(a)-(c)  with  G = L (12)  , * = t *,  and  n = LI  , and (2.1)  we  easily 


"*-V°  < C dfsup^{lld|l0  ||*-Lh*ll0  + II*  - uh,lo  llyd  - whydll0 


+ llgllo  IUa  - Wo}/lldllo 


From  (3.22)  we  get 


II  * - Ml  n - Ch  H tf'H  ' 


IUd  ' Wo^  11  V 2 


Using  (2.15)  and  (3.17)  we  have 


11  * - *h"o  - "h&»0  - cll)4  ■ V4"o  11  - Va"o 


ich  11*11,11^11, 
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Noting  that 


IU-VJI  • Ch  II  0(1  . 
h 0 — • 4 


11  dll ^ and  II y II 0 011  ^11  ^ , using  (3.  30)  and  combining  t 

obt  a in 
(3.  33) 
fox*  k * 1. 

Next  wo  consider  the  case  k ' 2.  Using  Theorems  3(a)  -(c)  with  G i.  (>;) 

and  n * I,  \ and  (2.i>  wo  have 
h d 

H V-livll,  i c *U1’  (lldll  II  t|- !'  (I’ll  ♦ llu-u  II  II  y - n y II  * |b(u  - n u,\  - I ’ 

h 0 a,  !.,<>:)  0 h 0 K Vh  0 '■<>  hW  o ' hV'  -i  h 

From  (t.22)  with  C « s*  * min (r-1 #k+ 1)  wo  got 

II*  ■ Ch*  II  Wll  - . 


"sinu  (2.1M  and  (1.17)  with  ( - s-l  wo  poo  that 

H\i  - W.K.  - Oh"  1 IU4I-  , 

' 'h  0 s+ 1 

and  usin.j  (1.17)  with  C - 1 and  (1.10)  wo  obtain 

II V , ■ ”,  y ,11  , ' Ch  If  dll  . 

\d  h.yd  0 — 0 

Finally  from  ( .1 . 7'0  with  v * y,  *•  - \ ^ , ( - 5-1,  and  ( 

lb(v  - V'd  - Vd’1  • C,'S  IUII5,1  lldll0  • 

Comtun \ng  those  or.t  \wM  oa  wo  have 


and  ( 3 . U ) wo 


( 3.  U> 


II  II  - (J-.  II  , ■ Ch  II  Ull  - , , 

h 0 s«l 


where  s - min(r-l,k*l)  and  k • 7 . 

Not  o that  (1.14)  a i von  an  improvement  ovoi  (i.i|)  only  foi  ki  1 , \.  r • 
(i.iM  improve  estimates  in  (71.  BabuSka  , Osborn,  and  V it  kar  ant  a ( M have  ptov. 
(1.  11)-(7.  U)  . 


d'  Kav  iat  t -Thtimas  mot  hod 

In  out  final  example  wo  study  a mixed  method  for  sooond  otdot  olliptio  probl 
introduced  by  Kaviatt  and  Thomas  1701.  Pot  <j  . I,,(il),  0 a convex  poison  in 

cons idot  the  modi’ 1 problem 

I 

-Ail  « u in  il 

K1  -1  0 on  1' . 

t ot  (l(diViil)  iy  , It. ,(11)1“:  div  V , i , (si)  ) with  t lie  norm 


i ' is' 


/ 

\ 


II  y II  , 

' (t(divtil)  ' 0 


- (Ilvli;  » lldiv  vir')'  •’ 

'0 


oims  w«' 


;j)  ' Midi! 


.i  t'  t hat 


(mat  «'• 

>d  ( t \) 

om : • 

IK*  , wo 


The  mixed  metluHt  wo  study  is  based  on  the  fol lowing  variational  formula! ion  of  (i.iM. 


KilUl  |l(dlV;s')  ' l Jl.')  Hill'll  tll.1t 

i'.  Ui>  y * y dx  ♦ v div  y dx  0 V v « H(iUvn>) 

i)  ii 

and 

(1.  I’M  J y*  (div  y « g)d*  0 V v*  < l . 

ll 

In  Theot  ero  l of  | .'01  it  is  shown  that  problem  ('.  U>)  -(  1.  17)  has  a unique  solid  ion 

( y . 4' ) II  bt  iv:  s’)  " l that  \t  is  the  solvit  ion  of  problem  1 h v ^ , ami  it  17i.nl  v. 

\n  add  it  ion  t lie  following  to. jnl.it  tty  insult  is  known  for  this  problem: 

If  •■)  < thon  g . iris,')  1 .111.1 

- 0 

( 1.  W>  IM  , oil  g II  % v g « 1 . («)  . 

One  i*.is  1 l y sees  that  tl.  *<0  - ( t . 17)  1 s an  example  of  problem  P with  V - ll(div;P), 

W - , \{  11  , til)  1 , a(y,vl  1 \i  • v ilx  , amt  b(y,v)  ~ f ij  div  u .lx 

h'  \) 

The  subclass  p of  data  tot  wh.ioh  (III)  is  satisfied  is  given  by  p - o s W * . siiuv 

a is  symmetric,  the  adjoint  problem  (2.*>),  (2,6)  with  0 - W - !,(»))  is  the  same  as 

probe Im  l*  and  thus  is  uniquely  solvable  fot  all  d . W*  . Iloneo  01.')  is  satisfied. 

Using  t t.  ) we  also  see  that  V . it  (ill  \ u (y\  , y - grad  \ , and 

d 0 ■tit  d 

(1.  V0  II  y .11 1 » l|\JI  i'll  dll  . 

M 1 d 0 


Wo  now  dosi't  t b*'  the  t mite  itimensional  sitbspaoos  used  m t hi'  approx imat  ion  sohomo . 
Pol  lowing  l-'Ol  we  bog  in  by  int t iMuo ing  the  space  ^ associated  with  the  unit  tight  tii 
angle  T in  the  (£.,  n)  -plane  whose  vettioes  ate  a ^ - (1,0),  a>  (0,1),  ,i;  (0,0). 

fot  k ' 0 an  even  intogoi  , def  ine  ^ 1 o be  the  spare  of  all  f vmot  ions  q of  the  form 

qj  - pol^(f. ,n)  ♦ n ♦ ...  ♦ a . C n , 


, . K»  1 . k A k/2*  l 

pol.(.‘. , n)  ♦ i'  u * i\.\n  . t.  n 

k 0 1 k . 


with 


k 

V (-1)  * (a.  t*  1 0 , 

i-0  1 

whet  e pol^t.  , *«)  denote-,  any  polynomial  of  degree  k in  the  two  vat  table1  • , n. 


For  K - l an  vxid  integer,  define  ^ to  be  the  space  of  all  functions  q of  the  form 


k+ 1 k 

*i1  • polk(i.,n)  » >>0f.  + »jf.  n 


. (Ml)/2  (k*l)/2 

(MU' 2 ' 


, . , k*l  „ , k , (k*  l)/2  (k+l)/2 

*2  ‘ P°lk(t'n)  * V'  * *ltn  4 •••  4 ^(k+ 1)  /2  L 

with 

(k*l)/2  (k+l)/2 

y (-D  v - i (-D ‘p  - o . 

i-0  i-0 

Now  consider  any  triangle  T in  the  (x^ ,x >) -plane  whose  vertices  are  denoted  by 

. * * 2 2 

a.  , l < i ' 3.  Let  V'  : x * F (x)  - B x b , B < Jf  ( m ) , b t IB  be  the  uniqu 
i — — T T TT  T T 

invertible  affine  mapping  such  that  F^la.)  - a ^ , 1 ^ i ' l . With  each  vector -valued 

function  v - (v^#v,)  defined  on  T , we  associate  the  function  y defined  on  T by 

1 -1 
Y ' J bt  Y • kt 

where 

dT  - det  (Bt)  . 

For  0 * h • l , assume'  that  \ is  a triangulation  of  made  up  of  triangles  T 

h 

whose  diameters  are  less  than  or  equal  to  h which  satisfy  the  minimal  angle  condition 


(see  Subsect  ion  La).  We  finally  consider  problem  with 


vh  * *Yt,  1 «(div»‘a  ! 


'h  ’ ”h 


• 0 ) 
T X.T 


" 1 v * H (divjT)  : v » £ ) 


W4  - U\  < L.  (tt)  : V T r 1 . v\  I « \\  ) , 

h h h h |T  k 

To  apply  out  t heot  ems  we  must  check  that  t hi'  appropriate  hypotheses  are  sal  isfied. 

Now  (HI)  is  trivially  satisfied  with  a • l.  In  the  proof  of  Theorem  I of  1201  it  is 

essentially  shown  that  there  is  an  operator  x : (H*(iJ))*  * satisfying 

1 2 

b(v  - n v , v' ) - 0 V v f III  (il)l  and  V *•  « W, 
h ' h 


Flirt  hetmore , tot  v . IH  (j.'l*.  i N 2 . we  have 


(3.40)  II v - nhvll0  <_  Ch*  llvll  ^ 1 < e ^min(r-l,  k+1), 

and 

(3.41)  II  div  (v  - tt  v ) II  < Ch™  lldiv  vll  , 0 < m < min  (r-2  , k+  1 ) . 

*v  h'v  0 — ^ m 

Using  the  regularity  results  (3.18)  and  (3.39)  we  easily  see  that  V [!l'l..)l~  ■<> 

that  (H5)  is  valid. 

We  next  observe  that  for  ^ t , div  ^h|T  < P^-  Hence  v^  < Z^  easily  implies 

div  v = 0 and  so  v.  t Z.  Thus  Z,  c z and  so  we  are  in  the  special  cases  of  Theorem: 
xh  Mi  h 

2 and  3.  Again  by  the  Corollary  to  Theorem  2,  (H4)  need  not  be  satisfied  in  order  to 
apply  (2.15),  so  we  shall  not  require  to  be  quasi-uniform. 

We  are  now  ready  to  derive  the  error  estimates.  Assume  that  41  t Hr((l)  , r ^ 3.  Prom 
(2.15)  and  (3.40)  we  obtain  for  k >_  0 , 

(3.42)  II u - u II  < Cllu  - null 

lh  0 - t h'v  0 


<_  ch*"  llsjlit  < ch1  Mt+1  , 


where  t = min (r-1 , k+1 ) . 

Now  applying  Theorem  3 (a) -(b)  we  get 

(3.43)  lli|>-i|i  H = sup  (b(y  - 

h ° deleft)  W 


hKd' 


-*) 


+ a(>ih  - - y + b(*  - v.'  vd  -n)  1/11  dllo 

for  all  y!,n  1 . Using  (3.41)  and  standard  approximat ility  properties  of  K , 

we  have 


(3.44) 


inf | b(y 

n«w  "d 


w.y  - > 4>-«P>  | < II  div  (y 
h.td  — jd 


y ,)  II  „ inf  Hwll 

>Kd  0 u 0 


< C lldiv  y II  hp  II i|>ll  , 
^d  0 u 


where  u = min(r,k+l),  and  clioosing  m = w - 2 in  (3.32), 


(3.45)  inf  |b(u  - n u,X.  - g)  I < II div (u  - n.  u) II  inf  II  X.  - nil 

1 n,  h*v  d 1 — % h'v  0 ..  d 


n.W, 


I'W,. 


< Ch1'-2  lldiv  ull  „ h2  II  x.ll 
— ■v  u-2  d 


provided  k - 1. 
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Us inv)  (2.1),  (2. IS)  . And  (3.40)  (with  < » u-l)  we  obtain  for  k 1 , 

u.4*>  !•%,  - V/  -hvd  - y I i - st«0  "\yd  - yd"0 

i c"v  - V'o  "\Xa  - Xa"o 

< rh1'-1  Hull  , h h y II  . 

- "c  U-l  ,'d  1 

Now  from  (3.43),  (3.44),  (3.4S),  (3.4b)  and  the  regularity  result  (3.39)  wo  obtain  ft 


(3.47)  ll*-i-  II  % v hM  11*11  , 

hO-  u 

where  \i  ■ min  (r.  Ml). 

To  obtain  an  estimate  when  k » 0,  we  choose  m - 0 in  (3.41)  to  obtain 

(3.48)  inf  |b(y  - - n)  | ^Hdiv(y  - inf  II  X ^ - nil  ^ 

n«w  1 n«w.  ‘ ( 

n h 

' c lldiv  yll0  till  \dll  j 

and  choose  i - l in  (3.40)  to  obtain  in  the  same  manner  as  in  (4.3i>)  that 

n.4-M  I«()ih  - v Va  - *d’  I i r,lV  - ”hV"o  11  \Xa  ‘ fc'o 

1rh  iiviiv  h II  . 

Cisnbining  ( 3.4  3),  ( 3.44)  with  k - 0 , (3.48),  ( 3.49),  and  the  requl.uity  renul*  (3.3‘M 


(3.  SO)  II  *-*^H  ) v Uhll  *11  > , k ■ 0 . 

We  note  that  estimate  (3.42)  was  obtained  in  (23,  IX  - 3.22al  and  that  (3.47) 
gives  an  improvement  over  the  result  in  l.M,  ix  - 3.22.11  in  the  cane  where  j h*  (el, 
y » ll r ( i23  , and  2 • r * Ml  . 
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